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We introduce a stochastic equation for the microscopic motion of a tagged particle in the single 
file model. This equation provides a compact representation of several of the system's properties 
such as Fluctuation-Dissipation and Linear Response relations, achieved by means of a diffusion 
noise approach. Most important, the proposed Langevin Equation reproduces quantitatively the 
three temporal regimes and the corresponding time scales: ballistic, diffusive and subdiffusive. 



I. INTRODUCTION 



in the following suggestive form 0] 



Since its introduction in 1965 due to Harris' pioneer- 
ing work [l| , the single file model (SF) has attracted more 
and more interest among the scientific community. Intro- 
duced first in the mathematical physics literature as an 
interesting though somewhat exotic topic, it has inspired 
over the last 40 years a large body of profound theo- 
retical studies j2] and detailed numerical investigations, 
including extensive Monte-Carlo and molecular dynamics 
simulations [3|. 

The motivations for the longstanding interest in this 
topic reside, on one hand, in its analytical tractability 
and, on the other hand, in its effectiveness as a descrip- 
tion of diffusion phenomena in real quasi one dimensional 
systems. As a matter of fact, since the direct observa- 
tion and manipulation of nanoscopic systems has expo- 
nentially evolved in the last decade, models suitable to 
account for the single particle diffusional mechanisms in 
constrained flow geometries, have been subjects of in- 
creasing attention. Remarkably, among these, the SF 
holds a preeminent position, since it correctly reproduces 
transport properties in a large category of quasi one di- 
mensional systems, where each particle is free to diffuse 
against its neighbors but is forbidden to overcome them 
1^. Transport processes of this type may be observed 
in nanoporous materials Q, in collective motion of ions 
through biological channels and membranes as well as 
in nanodevices and cellular flows [7] 

Along a mathematical line, the SF is perhaps the sim- 
plest interacting one dimensional gas one can consider: 
it consists of N unit-mass particles constrained to move 
along a line following a given dynamics. As the particles' 
interaction is purely hard-core, no mutual exchanges of 
the diffusants are allowed, i.e. they retain their order- 
ing over time {single filing condition). In spite of the 
intricacies of its mathematical derivation, the long time 
behavior of the single file dispersion relation can be cast 
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where p = j; (N — > oo, L —i- oo) is the file's density and 
(|A"|) is the absolute displacement of a non-interacting 
particle. If the free particle dynamics is characterized by 
a diffusivity D = ksT/^, then the relation ^ takes the 
form i 
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The predicted subdiffusive behavior has been reproduced 
experimentally in colloidal particles systems [3] and ob- 
served in molecular sieves (zeolites) [IQ]. We remark 
that the subdiffusive behavior represented by Eq. ([2]) 
will eventually be replaced by regular diffusion if for in- 
stance the particles are allowed to overtake each other 
[ll[, if there is only a finite number of particles [l2], or 
the particles move in a ring [Tsj . 

In Ref. lj| it was pointed out that the subdiffusive 
regime of a SF tagged particle occurs on the score of 
long ranged anticorrelations of its velocity, and/or of the 
jump's statistics of the collisional mechanism underly- 
ing its dynamics. Beside these persistent memory ef- 
fects, different mathematical derivations agree with the 
fact that asymptotically the tagged particle's probability 
distribution must be Gaussian with a variance growing in 
time according to ([2]) [2]. Together, these properties con- 
trast with the Continuous Time Random Walk (CTRW) 
scheme and its corresponding Fokker Planck representa- 
tion, for which a stretched Gaussian solution has to be 
expected (see RehfTs*] and references therein). Further- 
more, several subdiffusive systems in nature share the 
property of Gaussianicity with the SF, e.g. a monomer 
in a one dimensional phantom polymer [161] , the "translo- 
cation coordinate" of a two dimensional Rouse chain 
through a hole [17], a tagged monomer in an Edward- 
Wilkinson chain [ISj], de Gennes' defects along a polymer 
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during its reptation [l9| and solitons in the sine-Gordon 
chain 

In this paper we address the question of the micro- 
scopic effective description of the stochastic, anomalous 
motion of the tagged particle. Our aim is to extend the 
valuable Langevin approach, valid for the case of a dif- 
fusive Brownian walker, to the subdiffusive dynamics of 
a SF particle. We anticipate here that the Generalized 
Langevin Equation (GLE) [2l|,[23| provides the ideal the- 
oretical tool for such a goal, incorporating all the statisti- 
cal properties enjoyed by the particle. Within this frame- 
work the non-Markovian memory effects are achieved by 
means of a power-law damping kernel (23j , which is sim- 
ply added algebraically to the instantaneous friction of 
the surroundings. We note here that the GLE has re- 
cently been successfully used to describe several physical 
phenomena and market flows [2^ . 

The article is organized as follows: in sec|ll]we study 
the density profile dynamics by means of a Diffusion- 
noise approach and we connect the file's density fluc- 
tuations to the motion of a tagged particle. In seclIIII 
we introduce the GLE and we show the accuracy of the 
Langevin description by means of extensive molecular dy- 
namic simulations. 



II. DIFFUSION NOISE APPROACH 

We start by considering the file density dynamics. As 
stated inlH the system is composed of N Brownian point- 
like particles, all of unit mass, moving along a ring of 
length L and performing a stochastic motion according 
to the Langevin Equation (LE): 



Xi{t) = Vi(t) 



(3) 



where i G [^,N] denotes the particle's index, and the 
damping coefficient 7 and the random noise source ^i{t) 
satisfy the well-known Fluctuation-Dissipation Relation 
(often called second Kubo Theorem t22|]) (Cj(0>^j (^0) = 
2kBTjSijS{t — t'). The single filing condition turns 
out to be merely the interchange of two particles labels, 
whenever these suffer an (elastic) collision. However, as 
pointed out in [2^, all system properties which do not 
depend on particle labeling remain unchanged from those 
of an ideal gas, i.e. of N independent Brownian walkers. 

Let's first define the file density at a point x of the line 
at time t as 



n{x, t) 
dx 



(4) 



where n{x,t) refers to the number of particles in the bin 
[x — ^,x + ^] at time t. It is straightforward to note 
that the quantity p{x, t) is a local property of the file, 
independent of the relabeling of the particles due to col- 
lisions. A direct consequence of this is that the time 
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FIG. 1; (Color online). Panel (a): typical SF dynamics. 
The tagged particle (red) is subjected to the external peri- 
odic force F{t') = Fq cosut, while the other ones perform the 
usual stochastic dynamics. Panel (b): Linear Response re- 
lation (|18l) . Average velocity {v{t)) vs uit for fesT" = 1.0, 
7 = 0.5, p — 0.25 and different values of Fq and u. Three 
typical numerical curves (black lines) are fitted through the 
relation (I18|) . showing a linear dependence on the amplitude 
of the applied sinusoidal force. Averages have been taken 
over more than 2000 periods of each realization, performing 
at least 10 realizations corresponding to each uj. 



evolution of the file profile density can be described by 
the Diffusion noise equation for a 1 dimensional gas of N 
non-interacting Brownian particles [26j 



(5) 



having a recourse to the definition of a stochastic flux 
J{x,t) = —D-^p{x,t) +r]{x,t). The noise term ri[x,t) 
can be shown to be Gaussian and satisfy the following 
properties [27j : 



77(L,<) = 77(0,<) 
(77(x,t)) =0, 
{'q{x,t)ri{x' ,t')) = 2DS{x 



x')S{t^t'){p{x,t)), 



(6) 



the first one of which refers to the conservation of the 
particle's number {conserved noise) along the segment 
[0, L] with periodic boundary conditions. The remaining 
properties in ^ are required to fulfill the equations for 
the first two moments of p{x, t) 
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^{p{x,t))=D^{p{x,t)) 
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{p[xi,t)p(x2,t))c 



dt dx\ dx\ 
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where we made use of the short nota- 
tion {p{xi,t)p{x2,t))c = {p{xi,t)p{X2,t)) - 

{p{xi,t)){p{x2,t)). As it is apparent from the 
correlation function of the noise depends upon the 
particular solution of the diffusion equation in ([7]): in 
the following, as well as in the numerical simulations per- 
formed, we consider the case of a uniformly distributed 
file, namely (p(x,0)) = {p(x,t)) = p. 

The connection between the dynamics of the density 
over [0, L] and the the motion of a tagged particle in the 
single file system is achieved in the following way. Given 
two particle trajectories, xi{t) and X2(i), the number of 
particles between these has to remain constant in time 
because of the non-overlapping condition, this implies 



dt 



X2{t) 



p{x, t) dx = 0. 



(8) 



xi(t) 



Performing the derivative and making use of (O, the pre- 
vious relation then reads 

V2{t)p{x2{t),t) - J {X2{t),t) - Vl{t)p{xi{t),t) 

+Jixi{t),t) ^0; 

now, both the terms must set to zero irrespective of the 
particle labels: we can thus write down the equation for 
the single file particle as 



v{t)p{x{t),t) = J{x{t),t)^ 



(9) 



We notice that, although the relation ^ just introduced 
is exact, it is highly nonlinear. We are thus compelled 
to fall back upon two approximations in order to solve 
it. The first approximation is to assume that the density 
surrounding the particle position is essentially constant: 
p(0,t) = (p(0,t)) = p. The second approximation we 
put forward is to assume that the particles movements 
are correlated over a range equal to the displacement of 
the tagged particle, such that we can take the current 
J{x(t),t) to be equal to the current at the particles ini- 
tial position. Taking this to be at a: = we thus have 
J{x{t),t) = J(0, t). Notice that similar assumptions have 
been superimposed by Alexander and Pincus in a previ- 
ous treatment of the single file subdiffusive dynamics on 
a lattice [28}. 

With these approximations Eq.® gets the form 




FIG. 2: (Color online). Real (a) imaginary part (b) and phase 
(c) of p{lu) obtained through the relation (|18ll . for different 
Fo and lu. Simulation parameters are the same as in Fig[T] 
The dashed lines represent the formula (|15|) and are seen to fit 
quite well the asymptotic subdiffusive behavior of the mobility 
(see the inset of panel (b) where the low-frequency behavior 
of Qp{u!) is blown up). The initial diffusive behavior (orange 
lines) is responsible for the high-frequency regime uj ':$> ujd 
(dotted line in (a)). The GLE expression for the mobility 
(|45p furnishes a very good description of both regimes: solid 
black lines. The non-perfect agreement between theory and 
data is due to the small value of the ratio rd/rj, = 4 (see text 
and TABLE I). 



Moreover defining the Fourier transform (and its inverse) 
in the space and time domains as 



f{q,u;) = f+^dxdtf{x,t)e-^(^-~'^') 
the equation ([TU)) can be casted as 



(11) 
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J(q,w), 



which, by means of ^ and of the definition of the 
stochastic flux, reads 
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Let us now define the following quantity 
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using (|12p and the noise properties in ([5]), it is readily 
verified that the following equality holds: 



{v{u})v{uj')) =2kBTU[fi{uj)] 2n5{uj + uj') 
Furthermore a direct calculation of (fT3|l gives 



(14) 



(15) 



which, substituted in p4)l and consistent with the 
Wiener-Khintchine theorem, yields the asymptotic form 
of the velocity autocorrelation function (VAF) of a single 
file particle 



{v{t)vm = 



IknT 1 1 
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In passing from eq. lM]) to (fT6|) we imphcitly adopted 
the convention {v{t)v{0)) = {v{—t)v{0)), which uniquely 
determines the mobility /i(a') to be the Fourier-Laplace 
transform of the VAF, i.e. 



/x(w) 



knT 
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The relation (fT7| is known as the first Fluctuation- 
Dissipation theorem or Green-Kubo relation [23] ■ How- 
ever for linear systems the connection between corre- 
lation functions and mobility hinges on the domain of 
transport processes pof : within this framework the equa- 
tions and ([T5)) on one hand, and their statistical 
counterpart p7p on the other, encouraged us to explore 
the validity of a linear response relation of the type 



j{t)) ^ n Foe'' 



(18) 



Indeed, applying a periodic external force F{t) — 
Fq cos cut to a tagged particle and leaving the remain- 
ing surrounding ones unaffected (see Fig[lja)), the re- 
sponse (jl8p has been directly measured by meaning of 
extensive numerical simulations. The results are shown 
in Figlljb) and Fig[21 where the real, the imaginary part 
and the relative phase (p^j of the quantity /i(w) are dis- 
played. At first we note as the average velocity depends 
linearly upon the force amplitude Fq according to the 
equation (|18p . Most important, the low frequency be- 
havior of fully agrees with the analytical formula 
for the mobility fJ,{uj) given in (fT5|) (dashed line). We can 



interpret this result recalling that a particle undergoes a 
normal diffusive behavior up to a time scale = tj^j 
which is the time needed by a couple of particle to collide 
with one of the neighboring diffusants [l3 |: correspond- 
ingly we can assume that the tagged driven particle will 
feel the presence of the surrounding ones over frequencies 
smaller than the threshold Wd = ^ (dotted line). Con- 
versely, for CO ^ LOd, /i(w) coincide with the mobility of a 
free Brownian walker: 1/(7 — iuj). The numerical results 
in Fig[2]can thus be summarized by writing, besides the 
relation (fT8|) . the mobility fi{uj) as 




(19) 



The numerical evidence of the effectiveness of a linear 
response relation allows us to rewrite the eg. p^ as 



(20) 



where the introduced noise satisfies {^{t)) ~ 0. On the 
other hand, its spectrum 5^(w) exhibits two different 
regimes according to (111 



2kBT-i 



uj :^ u!d 

(21) 

The file particles surrounding the tagged one thus act as 
an additional bath responsible for the onset of subdiffu- 
sional behavior. The nature of this long-ranged corre- 
lations in the noise source can be easily understood in 
terms of the collisional interaction between the file com- 
ponents. Note in fact that the expression (I^T]) leads to a 
slowly decaying positive correlated noise {S,{t)£,{Q)) oc -^: 
this non-trivial finding is the signature of constrained ge- 
ometry systems. Indeed, although the collisions tend to 
tie back the particle motion, leading to the negative ve- 
locity correlations (|16p. the noise provides the way to 
maintain {v{t)) = on time scale of the order I/7. An- 
other way to say this, is that in a collision a particle 
exchanges velocity and noise. Anyway, this is a manifes- 
tation of the Fluctuation-Dissipation theorem. Further- 
more, the non-Markovian power-law nature of the noise 
spectrum characterizes the asymptotic Fractional Brow- 
nian Motion (FBM) of the tagged particle [s?]. 

We end this section stressing that, using the formalism 
so far developed, we can calculate the other asymptotic 
statistical properties of the system. In fact it is straight- 
forward to write down an expression for the particle's 
position similarly to what we did in (jl2p for its velocity: 
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For instance, making use of and (fT^ we get 



(22) 



5 



150 



{x{u;)v{u;')) 



27 



: 2tt5 (w ■ 



(23) 



III. GENERALIZED LANGEVIN DESCRIPTION 

In this section we will collect the results outlined in the 
previous one and will put them in a consistent compact 
formulation. We emphasize three fundamental properties 
that such a representation must incorporate 

i) the Linear Response relation (LR) must hold: ([18]); 

ii) the Fluctuation Dissipation theorem (FDT) is also 
vahd: dUD; 

in) all the statistical properties exhibit two different 
behaviors, Brownian or subdiffusive motion, depending 
on whether time is smaller or larger than r^. We will 
in the following assume that the time scale r;, = I/7 is 
smaller than such that the particles are moving diffu- 
sively (in contrast with ballistically) before colliding with 
each other. 

Our aim is thus to write down an effective equation for 
the microscopic dynamics of a tagged particle in single 
file systems, according to i), ii), Hi). Such an equation 
turns out to be 
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FIG. 3: (Color online) Second moments of observables: 
Rescaled VAF {v(t)v(0)) (a) and position-velocity correlation 
function {\x{t) — x(0)\ v{t)) (b) vs t/rd for several densities p 
and damping 7. ksT = 1.0. Dashed lines are the asymptotic 
formulae (I16|l and 1)23^ respectively, whereas the solid ones 
represent the CLE predictions (|35|l and (|40|) for p = 0.05 and 
7 = 5.0. To improve our statistics we averaged over 3000 
different realizations. 
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Several new symbols have been introduced in the pre- 
vious expression deserving an explanation. Firstly the 



The noise appearing in (|24p satisfies the properties 
(I^Tj) . Indeed we recall that the diffusive microscopic time 
scale Til = - can be expressed as t;, = ^ so that the 



plays the same role as Td in the equation ^ can be re-casted as 



quantity Ud = / 
collisional representation of the particle's motion [lj|: it 
accounts for the number of collisions a particle suffers be- 
fore attaining subdiffusive behavior. Secondly the quan- 
tity 7i/2: which has the dimension of [l/t^/^] , is the gen- 

eralized damping coefficient and is equal to ( — 1 since 

Td is the unique relevant time scale of the system. The 
third symbol in ()24ll is the Riemann-Liouville fractional 
operator (soj : 



v{t) 



j{t-t')v{t')dt' + i{t), 



(27) 



introducing the generalized damping and defining the S 
function to contribute half at the end point of an integral: 
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Thence such a definition allows us to express the proper- 
ties in ((2T|) through the compact notation of the Gener- 
alized Fluctuation- Dissipation theorem [22i\ 



Note that if we use the definition of the Caputo fractional 
derivative [Sll 
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then ([24|) takes the form 

d^l'^x{t) 



dt\ 



(29) 



The solution of the equation can be easily achieved 
by means of the Laplace transform: 



d'^x{t) ^^dxlt) 



x{s)='^+v{Q)i;{s)+~i{s)i^{s) 
v{s) ^ v{0)fL{s) + ~^{s)fL{s), 



(30) 



dt 
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where the mobility jl{t) in the domain is given by 
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and 
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It is immediate to verify that the expression in (I3ip 
matches the predicted two regimes in p^ . In the time 
domain the solution of the Generalized Langevin Equa- 
tion in ([M)) takes the simple form 



x{t) = x{0) + viO)i>it) + /o |(t')^^(i - t') dt' 
vit)^v{0)j,{t) + j*iif)fl{t-t')dt', 



(33) 



from which it is apparent that the joint probability dis- 
tribution for X and u is a Gaussian in this formulation. 
Position and velocity are in fact linear functionals of 
^(t) which is a (non-Mar kovian) Gaussian random pro- 
cess. Such a property corroborates all the previous an- 
alytical derivations for the probability distribution of a 
tagged particle moving subdiffusively in a single file sys- 
tem (see for example Ref.[2] and references therein) and 
it is clearly at odds with a Fractional Fokker-Planck de- 
scription of the corresponding subdiffusive process, which 
leads to a stretched Gaussian solution T^. 
Furthermore we point out that the GLE ([M)) correctly de- 
scribes the time behavior of all the observable moments 
of position and velocity, not only in the long time asymp- 
totic regime but also at the initial diffusive stage. To see 
this it is sufficient to start from the general solutions ([55]) 
and put in some physical assumptions [s^ . For instance 
the first moments of velocity and position will read 



(xit)) = {xm + {vmm 
{v{t)) = {vmr^t). 



(34) 



The second moments, however, are more interesting 
quantities: for the VAF, assuming (i;^(0)) — ksT, it is 
straightforward to prove 



C(t) = (^^(^^(O)) - keTfiit) 



(35) 



which is the Generalized First Kubo theorem. The nu- 
merical evidence of this is given in Fig|31[a) , where several 
rescaled curves are plotted against the analytic function 
in (j35p after numerical inversion of the mobility in (|3ip . 

The excellent agreement between the analytical de- 
scription yielded by (|24p and the numerical data is even 
more apparent looking at the mean square displacement 
of the tagged particle (Fig[3|). Indeed from ([55)1 one ob- 
tains 
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FIG. 4: (Color online) Mean square displacement 5a;^(f) = 
\xi{t) — a;i(0)]'^ I'M for different values of p and 7 with 
ksT — 1.0 plotted versus t/rd. Data have been rescaled by 
x\ = Ij (p^\/vr) y ^^is [1j1- The black solid lines repre- 
sent the theoretical predictions as given by the GLE formula 
(|36p : they provide an excellent description of the three dif- 
fusive regimes. Data have been averaged over 30 different 
realizations for A'^ of the order of 10^. 



8x^(t) = l\x(t) - x(Ssf) = 2 {v^{0)) / ^A(i') dt', (36) 

Jo 

where it is possible to recognize the theorem stated in 
Ref. [11] provided that 



(37) 



The exact expression given in (|36p quantitatively repro- 
duces the three stages of the Sx^ (t) curves in FigH] (bal- 
listic, diffusive, subdiffusive) as well as the characteristic 
time scales (r?, and t^) on which the crossovers between 
them take place 



Sx'^it) 



ksTt^ t<^Tb 
2Dt n^t^Td 

'^X^V't t > Td. 
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An analytical inversion of the formula in (I32p including 
the crossover to the subdiffusive regime can be achieved 
by neglecting the inertial (ballistic) term leading to 
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In general we believe the GLE (|24p provides a very good 
description of the tagged particle's stochastic motion only 
when the condition T}, ^ is fulfilled (see TABLE I). 
In other words the particle must attain a truly diffusive 
regime before getting a collision with one of its nearest 
neighbors. One can see that the proposed GLE cannot 
work well in the opposite case where Tf, > since all 
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TABLE I: Values of Tb, ra and their ratio for the values of p 
and 7 used in Figs. El O H and [5] 



interactions with neighboring particles vanish according 
to Eq. (p4|) as 7 ^ 0. This cannot be true, since the 
particle will still collide and exchange momentum with 
its neighbors even though the rest of the friction with 
the surroundings vanish. We have not found a GLE that 
works well in the case Tf, > t^. 

The close comparison between theory and numerics when 
Tb <C Td is also displayed by the position-velocity correla- 
tion function {x{t)v{t)) for which the following expression 
holds 



(40) 



In Fig[3lb) we compare the numerical data with both 
(j40| and ([23ll , which is expected to work well only in the 
asymptotic regime. 

The relaxation of the second moment of the velocity 
{v^{t)) to the asymptotic value ksT is given by 

{v'{t)) = {v'm f,'{t) + ksT [1 - ^,\t)] , (41) 

and deserves particular attention. In Ref. [1^ it has 
been shown that the Fractional Brownian process (3^ 
generated by a Fractional Langevin equation and the 
stochastic process corresponding to the relative Frac- 
tional Kramer equation (3^ are on average the same, 
except for the second moment of the velocity. We point 
out that such a discrepancy is a common problem in 
non-equilibrium statistical mechanics, whenever one is 
concerned to pass from a Generalized Langevin dynam- 
ical description of a non-Markovian process [21] to the 
corresponding Fokker-Planck equation for the probabil- 
ity distribution [ssj (see for instance the discussion in 
Ref. [131 on the inconsistency of a retarded Fokker-Planck 
equation of the Rubin model [1^). In the GLE ^ the 
{v'^{t)) relaxation process is dictated by the Brownian 
dynamics (see Fig. [5]) as long as n <^ Td- The Fractional 
Kramer equation corresponding to would read 



dP{x,v,t) _^^dP{x,v,t) 



dt 
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LFpP{,x,v,t), (42) 



where Lpp = -^v + kpT-^. In the Tf, <C Td limit the 
(ti^(i)) FKE expression yielded by is the one provided 




FIG. 5: (Color onhne). Second moment of the velocity {v^{t)) 
vs t, displayed for several densities p and damping 7. ksT = 
1.0. The initial tagged particle velocity «(0) has been set to 
0. The expression (|41|l (solid black lines) shows the typical 
Brownian exponential relaxation to fcsT: for all the displayed 
curves Tb/rd 1. To improve our statistics we averaged the 
numerical curves over 3000 realizations. 



by the usual Fokker-Planck, i.e. (w^(t)) fke — ksT e^^'^* 
once v(0) = 0, in agreement with the exponential satu- 
ration shown in Fig. O This could lead to believe that 
the GLE description (|24p and the Kramer Eq. ([1^ are de 
facto equivalent in the range Tf, <^ Td- Nevertheless we 
expect that the P{x, t) solution of (j42p is still a stretched 
exponential instead of a Gaussian: this fundamental dif- 
ference casts some general doubt on the possibility to 
determine a Fokker-Planck equation for systems whose 
microscopic dynamics is represented by a FLE. 

The last part of this section is devoted to the prop- 
erty i), namely the validity of LR. In the presence of an 
external force acting from time equal zero and onwards, 
(I24|-(l27l) takes the form 



i(t) = v{t) 

m = - Jo 7 (i - 1') v(t') dt' + F{t) + m. 



(43) 



The LR ([T^ in s— space can thus expressed by using ([31 
as 



{v{s)) ^ fi{s)F{s), 



(44) 



so that , thanks to the expression ipi]) , the real and imag- 
inary part of the mobility jlit) read 



3?[A(^)] = 



vi 



-1/2 



{'t+V2nl 7l/2a>-i/2)^ + (V2«^^i/2C^-i/2_^)^ 

cx,r~/ \-i _ i^-y^ng 71/2'^"'"^^ 

^ ^^^^'^ ^ (7+V2 n3 7l/2'^-i/^)'-H(\/2n5 j,/^uj-^/^-uY 

(45) 

Remarkably, in Fig. [5] both functions ([l5)) (solid black 
lines) are shown to fit quite well the outcome of our nu- 
merics. 

Finally in the case of a constant external force, i.e. 
F{t) = for t > in (gS]), the LR provides that the 
drift {x{t)) p satisfies the Generalized Einstein relation 
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(x(0)i. = (x(0))+F— ^ (46) 

that has been shown to be vahd on file systems moving 
on a ring [38,] and on an infinite one-dimensional lattice 

IV. CONCLUSION 

In this paper we introduced the effective equation rul- 
ing the microscopic stochastic motion of a SF tagged par- 
ticle. Starting from a diffusion-like formalism for the file 
density dynamics, we were led to several properties for 
which the GLE (l^^ - (P7)) can be regarded as a representa- 
tion. Indeed the GLE formalism provides an elegant rep- 
resentation of Generalized Fluctuation-Dissipation the- 
orems (|29p - ([55)) . Linear Response pi)) and Generalized 
Einstein relation (|46)) . all remarkable properties satisfied 
by the particle, both in its diffusive and in subdiffusive 
phase. 

Nevertheless we want to stress that the GLE is just a 
description, though very good, of the stochastic motion 
of the SF particle and, in this perspective, is valid within 
certain limits. However, along the same line, within cer- 
tain approximations the Langevin Equation provides an 
excellent effective equation for the motion of a Brown- 
ian particle immersed in a thermal bath, i.e. it describes 



well the ballistic and diffusive regimes of the particle, 
but one would not expect it to be exact in the region 
of crossover between the two regimes. In SF systems 
the tagged particle is subjected to two types of thermal 
baths, as is apparent from (j24[) : the first is mimicked by 
the usual Markovian uncorrelated noise, whereas the sec- 
ond, physically embodied by the surrounding file's par- 
ticles, is achieved by the introduction of an additional 
non-Markovian term, responsible for the strong memory 
effects. Surprisingly this "sum of thermal baths" turns 
out to be well described (in the same manner as for the 
Brownian particle/LE) by simply the algebraic sum of 
two independent terms in the stochastic equation of mo- 
tion. Alternatively one could view (f24|) as an usual LE 
([5]) where the relabeling-coUisional symmetry accounts 
for the fractional term. 
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